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Abstract 

We derive the elastic properties of a cylindrical cloak for in-plane 
coupled shear and pressure waves. The cloak is characterized by a 
rank 4 elasticity tensor with 16 spatially varying entries which are 
deduced from a geometric transform. Remarkably, the Navier equa- 
tions retain their form under this transform, which is generally untrue 
[Milton et al., New J. Phys. 8, 248 (2006)]. We numerically check 
that clamped and freely vibrating obstacles located inside the neutral 
region are cloaked disrespectful of the frequency and the polarization 
of an incoming elastic wave. 

Recently, significant progress has been made on the control of acoustic 
and electromagnetic waves. Transformation based solutions to the conductiv- 
ity and Maxwell's equations in curvilinear coordinate systems, subsequently 
reported by Greenleaf et al. [1] and then by Pendry et al. [2] and Leon- 
hardt [3], enable one to bend electromagnetic waves around arbitrarily sized 
and shaped solids. More precisely, the electromagnetic invisibility cloak is 
a metamaterial which maps a concealment region into a surrounding shell: 
as a result of the coordinate transformation the permittivity and permeabil- 
ity are strongly heterogeneous and anisotropic within the cloak, yet fulfilling 
impedance matching with the surrounding vacuum. The cloak thus neither 
scatter waves nor induces a shadow in the transmitted field. In [3], a cylin- 
drical electromagnetic cloak constructed using specially designed concentric 
arrays of split ring resonators, was shown to conceal a copper cylinder around 
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8.5 GHz. The effectiveness of the transformation based cloak was numeri- 
cally demonstrated solving Maxwell's equations using finite elements for an 
incident plane wave (far field limit) j5] and for electric line current and mag- 
netic loop sources (near field limit) [6]. In [7j, a reduced set of material 
parameters was introduced to relax the constraint on the permeability, nec- 
essarily leading to an impedance mismatch with vacuum which was shown 
to preserve the cloak effectiveness to a good extent. Other routes to invis- 
ibility include reduction of backscatter [SJ and cloaking through anomalous 
localized resonances, the latter one using negative refraction 0, [10]. To date, 
a plethora of research papers has been published in the fast growing field of 
transformation optics. 

However, transformation based invisibility cloaks applied to certain types 
of elastodynamic waves in structural mechanics received less attention, since 
the Navier equations do not usually retain their form under geometric changes 
[TIj [T2] . One simplification occurs for cylindrical geometries, whereby out- 
of-plane shear waves decouple from in-plane waves. However, in-plane shear 
and pressure waves remain inherently coupled. Earlier proposals for neutral 
inclusions include using asymptotic and computational methods to find suit- 
able material parameters for coated cylindrical inclusions [13] . The latter has 
proved successful in the elastostatic context in the case of anti-plane shear 
and in-plane coupled pressure and shear polarizations. However, neutrality 
breaks down for finite frequencies. 

Other avenues to elastic cloaking should therefore be investigated. For in- 
stance, Cummer and Schurig demonstrated that acoustic waves in a fluid un- 
dergo the same geometric transform as electromagnetic waves do and there- 
fore retain their form [14]. This result has been since then generalized to 
three-dimensional acoustic cloaks for pressure waves [TBI US] • Importantly, 
such cloaks require an anisotropic mass density which can be obtained via 
a homogenization approach, which presents the advantage to be broadband 
|17j . Acoustic cloaking for linear surface water waves was chiefly achieved 
via the same mechanism in between 10 and 15 Hertz [IB] . 

In the present letter, we show that it is also possible to design a cylindrical 
cloak for in-plane coupled pressure and shear elastic waves. We demonstrate 
theoretically its unique mechanism and further perform finite element com- 
putations checked again analytical calculations of the Green's function for 
the Navier equations in transformed coordinates. The main difference with 
previous work [9] is that our elasticity tensor in the transformed coordinates 
is no longer symmetric, which is a necessary condition for the Navier equa- 
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Figure 1: Geometric transform of eqn.([2]) from (r, 9) (Fig{T](a)) to (r',9') 
(Figfjjb)); tq and t\ are the inner and the outer radius of the cylindrical 
cloak, respectively. The elastic constitutive tensor and the density in the 
undeformed and in the deformed domains are denote by C, p and C', p', 
respectively. 



tions to retain their form. Quite remarkably, we find that the density remains 
a scalar quantity in the transformed coordinates. 

We consider the in-plane propagation of time-harmonic elastic waves gov- 
erned by the Navier equations 

V • C : Vu + pcAi + b = , (1) 

th 

where u is the displacement, p the density, C the 4 -order constitutive 
tensor of the linear elastic material and b = b(x) represents the spatial 
distribution of a simple harmonic body force b(x, t) = b(x) exp(iut), with uj 
the wave-frequency and t the time. 

We introduce the geometric transform (r, 9) — > (r', 9') of [U [2] 

J r' = r + r -^r , 9' = 9 , for r < rl 

\ r' — r , 9' = 9 , forr > rl 

shown in Fig. (OQ) and expressed in cylindrical coordinates r = \/ x\ + x\ and 
9 = 2atan(x2/(xi + \/ x\ + x\)), with ro and r\ the inner and outer radii of 
the circular cloak, respectively. 

By application of transformation ([2]), in the region r' G [ro, ri] the Navier 
equations ([I]) are mapped into the equations 

V • C : Vu + p'uj 2 m = , (3) 
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Figure 2: Elastic cloak in an elastic medium subjected to a concentrated load, 
(a) Displacement magnitude u = -\fu\ + w|; (b) deformation en = |^-; (c) 
deformation e 22 = ; (d) deformation e 12 = e 21 = ±(gi + fsa). 

where the body force is assumed to be zero, the stretched density p' = 
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) 2 p and elasticity tensor C' has non zero cylindrical components 



c: 




(4) 



with A and p, the Lame moduli characterizing the isotropic behavior described 



Interestingly, the transformation ([2]) preserves the isotropy of the density, 
which remains a scalar (yet spatially varying) quantity in (J3]), and avoids any 
coupling between stress and velocity. This is a very unlikely situation for elas- 
todynamic waves propagating in anisotropic heterogeneous media [11]. We 
also note that the proposed formulation poses no limitations on the applied u> 
ranging form low to high frequency, as the elasticity tensor does not depend 
upon uj. 

We report the finite element computations performed in the COMSOL 
multiphysics package. The elastic cloak of equation (j3J) is embedded in an 
isotropic elastic material with Lame moduli A = 2.3 and p, — 1 and density 
p — 1, which are realistic normalized parameters corresponding to fused 
silica. The elastic cloak has inner and outer radii r = 0.2 m and T\ = 0.4 m, 
respectively. 

The system is excited by an harmonic unit concentrated force applied in 
direction x\ and vibrating with angular frequency u = 40 Hz. A perfectly 
matched cylindrical layer has been implemented in order to model the infinite 
elastic medium surrounding the cloak (cf. outer ring on panels a, b, c and d 
of Fig f5]); this has been obtained by application of the geometric transform 



where r 2 = 1 m is the inner radius of the outer ring in Fig. ([2]). 

In Fig. ([2]) we clearly see that the wave patterns of the displacement and 
deformations are smoothly bent around the central region within the cloak 
(where the magnitudes are nearly zero). Although the coupling of shear and 
pressure waves generated by the concentrated force creates the optical illu- 
sion of interferences, the comparison with the harmonic Green's function in 
homogeneous elastic space (see [20] ) reported in Fig. (j3J) shows, at least qual- 
itatively, that there is neither forward nor backward scattering. The absence 
of scattering is better detailed in Fig. H] where results of Fig. (j2J) and Fig. (JSJ) 
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Figure 3: Harmonic Green's function in homogeneous elastic space, (a) Dis- 
placement magnitude u; (b) deformation e n ; (c) deformation e 2 2] (d) defor- 
mation £i 2 . 

are compared: the perfect agreement of the displacement fields in the exter- 
nal matrix with and without the cloak is shown, the distortion being bounded 
to the central region delimited by the cloak. These are non-intuitive results, 
as the profiles of the horizontal and vertical displacements in Fig. H] should 
display a visible phase shift, since the associated acoustic paths are different. 
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More precisely, let us look at the expression of the elasticity tensor given in 
(jlj). On the inner boundary of the cloak, that is for r = r , its components 
C' rrrr and C r0r9 vanish, whereas its components C 9gdg and C 9rdr tend to in- 
finity. This physically means that pressure and shear waves propagate with 
an infinite velocity in the 9 direction along the inner boundary, which results 
in a vanishing phase shift between a wave propagating in a homogeneous 
elastic space and another one propagating around the concealed region: this 
explains the superimposed profiles of horizontal and vertical displacements 
in Fig. HI 
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Figure 4: Comparison between numerical results in presence of the elastic 
cloak of Fig. [2J (black dots) and Green's function in homogeneous elastic 
space of Fig. [3] (grey lines). Results are given along the line AB detailed in 
Fig. [2(a). (a) Horizontal displacement u\\ (b) Vertical displacement u 2 . 

In conclusion, we have proposed an elastic cloak bending the trajectory 
of in-plane coupled shear and pressure elastic waves around a cylindrical 
obstacle. The cloak can be designed by the use of heterogeneous density 
and heterogeneous and anisotropic elastic stiffness; the distribution of the 
physical properties has been obtained with the introduction of stretched co- 
ordinates. Our results open new vistas in cloaking devices for elastodynamic 
waves in anisotropic media, yet with an isotropic density. 

Part of this work was performed while MB was working under the EPSRC 
grant EP/F027125/1. 
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